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Abstract 

The phenomenon of neutrino oscillations has been established as the leading mechanism behind 
neutrino flavor transitions, providing solid experimental evidence that neutrinos are massive and 
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I. INTRODUCTION 

Since the results of the Super-Kamiokande experiment in Japan in 1998 [lj, the phe- 
nomenon of neutrino oscillations has been established as the leading mechanism behind 
neutrino flavor transitions. This result was followed by a first boom of results from several 
international collaborations (e.g. SNO, KamLAND, K2K, MINOS, and MiniBooNE) on the 
various neutrino parameters. Certainly, these solid results have pinned down the values on 
the different parameters to an incredible precision given that neutrinos are very elusive par- 



ticles and the corresponding experiments are extraordinarily complex (see Ref. [2]). Never- 
theless, it is a fact that the present Standard Model (SM) of particle physics is not the whole 
story and needs to be revised in order to accommodate massive and mixed neutrinos, which 
leads to physics beyond the SM. With the upcoming results from the running or future neu- 
trino experiments (e.g. Daya Bay, Double Chooz, RENO, KATRIN, T2K, NOz/A, OPERA, 
ICARUS), there will be a second boom of results, and we will hopefully be able to determine 
the missing neutrino parameters such as the leptonic CP-violating phase (important for the 
matter-antimatter asymmetry in the Universe), the sign of the large mass-squared difference 
for neutrinos (important for the neutrino mass hierarchy), and the absolute neutrino mass 
scale, but also the next-to-leading order effects in neutrino flavor transitions. 

In future neutrino experiments (and in particular for a neutrino factory, /3-beams, or 
superbeams), new physics beyond the SM may appear in the form of unknown couplings 
involving neutrinos, which are usually referred to as non-standard neutrino interactions 
(NSIs). Compared to standard neutrino oscillations, NSIs could contribute to the oscillation 
probabilities and neutrino event rates as sub-leading effects, and may bring in very distinctive 
phenomena. Running and future neutrino experiments will provide us with more precision 
measurements on neutrino flavor transitions, and therefore, the window of searching for 
NSIs is open. In principle, NSIs could exist in the neutrino production, propagation, and 
detection processes, and the search for NSIs is complementary to the direct search for new 
physics conducted at the LHC. 

The concept of NSIs has been introduced in order to accommodate for sub-leading effects 
in neutrino flavor transitions. Previously, alternative scenarios for neutrino flavor transitions 
such as neutrino decoherence, neutrino decay, and NSIs have been studied (see below for 
references), but now, such alternatives are only allowed to provide sub-leading effects to 



neutrino oscillations. In the literature, there exist several theoretical and phenomenological 
studies of NSIs for atmospheric, accelerator, reactor, solar, and supernova neutrinos. In 
addition, some experimental collaborations have obtained bounds on NSIs. 

This review is organized as follows. In Section [Til we introduce the concept of neutrino 
flavor transitions with NSIs. First, we present standard neutrino oscillations, and then, 
we consider other scenarios for neutrino flavor transitions including NSIs. At the end of 
this section, we discuss so-called NSI Hamiltonian effects of neutrino oscillations. Then, 
in Section IIII[ we describe NSIs with three neutrino flavors, since there are at least three 
flavors in Nature. Especially, we consider production, propagation, and detection NSIs 
including the so-called zero-distance effect. In addition, we present mappings for NSIs and 
approximate formulas for two neutrino flavors that can be useful in some settings. Next, in 
Section IIV} we study different theoretical models for NSIs including e.g. see-saw models. In 
Section |V] we investigate the phenomenology of NSIs for different types of neutrinos such 
as atmospheric, accelerator, reactor, solar, and supernova neutrinos, whereas in Section IVI| 
we review phenomenological bounds on NSIs. In addition, in Section IVIH we examine 
experimental sensitivities and the future discovery reach of NSIs. Finally, in Section IVHH 
we present a summary and state our conclusions. 

II. NEUTRINO FLAVOR TRANSITIONS WITH NSIS 

In this section, we present the basic ingredients for neutrino oscillations (based on the two 
facts that neutrinos are massive and lepton flavors are mixed), which is the leptonic mixing 
matrix, a Schrodinger-like equation for the evolution of the neutrinos, and the values of 
the fundamental neutrino oscillations parameters, i.e., the neutrino mass-squared differences 
and the leptonic mixing parameters. Then, we discuss some historic alternative scenarios for 
neutrino flavor transitions. Finally, we study NSI Hamiltonian effects of neutrino oscillations. 

A. Neutrino oscillations 

Indeed, there are now strong evidences that neutrinos are massive and lepton flavors are 
mixed. In the SM, neutrinos are massless particle, and therefore, the SM must be extended 
by adding neutrino masses. The lepton flavor mixing is usually defined through the leptonic 
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mixing matrix U that can be written as 
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which relates the weak interaction eigenstates and the mass eigenstates through the leptonic 
mixing parameters (9 12 , #13, # 23 , 5 (the Dirac CP- violating phase), as well as p and a (the 
Majorana CP- violating phases). In the so-called standard parameterization, U is given by 
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where = cos(%) and s^- = sin(%) 



The time evolution of the neutrino vector of state z/ = \ p e u T j describing neutrino 
oscillations is given by a Schrodinger-like equation with a Hamiltonian H, namely, 



i^ = ^[MMt + diag(A0,0 



Hu. 



(3) 



where E is the neutrino energy, M = U diag(mi, m 2 , m 3 ) U T is the neutrino mass matrix, 
and A = 2y/2EGpN e is the effective matter potential induced by ordinary charged-current 
interactions with electrons 3|, |4|. Here, m\, m 2 , and m 3 are the definite masses of the neu- 
trino mass eigenstates and N e is the electron density of matter along the neutrino trajectory. 
Quantum mechanically, the transition probability amplitudes are given as overlaps of dif- 
ferent neutrino states, and finally, neutrino oscillation probabilities are defined as squared 
absolute values of the transition probability amplitudes. Thus, flavor transitions occur dur- 
ing the evolution of neutrinos. For example, in a two-flavor illustration (in vacuum) with 
electron and muon neutrinos, a neutrino state can be in a pure electron neutrino state at 
one time, whereas it can be in a pure muon neutrino state at another time. 

Using global fits of data from neutrino oscillation experiments (and especially, results 
from long-baseline neutrino experiments including matter effects), the values given in Table |T] 
have been obtained for the fundamental neutrino oscillation parameters [5|. Note that these 
values have been found without taking sub- leading effects such as NSIs into account. Open 
questions that still exist are: Are neutrinos Dirac or Majorana particles? What is the 
absolute neutrino mass scale? What is the sign of the large mass-squared difference Am^? 
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TABLE I. Present values of the fundamental neutrino oscillation parameters. See also Refs. [6j, [7J 
for two other global analyses. 

Is there leptonic CP violation? Do sterile neutrinos exist? However, recently, one has also 
been concerned with the following two questions in the literature: Are there NSIs? Is there 
non-unitarity in leptonic mixing? The intention of this review is to bring some insight into 
these last two questions (with emphasis on the first question). 



B. Other scenarios for neutrino flavor transitions 



The Super-Kamiokande, SNO, and KamLAND experiments have provided strong ev- 
idence for neutrino flavor transitions and that the theory of neutrino oscillations is the 
leading description [l, 8, 9]. In various neutrino oscillation experiments, precision measure- 
ments for some of the neutrino parameters, i.e., Am 2 1; [Am 2 ^, 612, #13, and 623, have been 
obtained, whereas other parameters are still completely unknown such as sign(Am| 1 ) and 
5, as well as the Majorana CP-violating phases and the absolute neutrino mass scale. 

Other mechanisms could be responsible for flavor transitions on a sub-leading level (see 
e.g. Ref. ID]). Therefore, we will phenomenologically study "new physics" effects due to 
NSIs. In the past, descriptions for transitions of neutrinos based on neutrino decoherence and 
neutrino decay have been extensively investigated in the literature lll-l53|. However, now, 



such descriptions are ruled out by available neutrino data as the leading-order mechanism 
behind neutrino flavor transitions 
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54j-l57|. but these descriptions could still provide sub- 



leading effects. In what follows, I will not consider neutrino decoherence and neutrino decay, 
but instead focus on NSIs, which are interactions between neutrinos and matter fermions 
(i.e., u, d, and e) that additionally affect neutrino oscillations, as a sub- leading mechanism 
for neutrino flavor transitions. 



5 



C. NSI Hamitonian effects of neutrino oscillations 



In general, NSIs can be considered to be effective additional contributions to the standard 



vacuum Hamiltonian H that describes the neutrino evolution (see e.g. Ref. 58]). Thus, any 
Hermitian non-standard Hamiltonian effect H' will alter the original Hamiltonian into an 
effective Hamiltonian: 

H eS = H + H'. (4) 
For example, neutrino oscillations in matter are described by 

H' = H mattCT = ^diag(A, 0, . . . , 0) - -±=G F N n t n , (5) 

where Gp is the Fermi weak coupling constant, N n is the nucleon number density, and t n is 
the n x n unit matrix. Just as the presence of matter affects the effective neutrino param- 
eters, the effective neutrino parameters will be affected by any non-standard Hamiltonian 
effect. The non-standard Hamiltonian effects can alter both the oscillation frequency and 
;he oscillation amplitude and they can be classified as "flavor effects" or "mass effects" 
Furthermore, non-standard Hamiltonian effects (such as NSIs) will lead to resonance 



conditions 
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which are modified (and more general) versions of the famous Mikheyev- 



Smirnow-Wolfenstein effect |3|,y,|59|. 



III. NSIS WITH THREE NEUTRINO FLAVORS 



The phenomenological consequences of NSIs have been investigated in great detail in the 



literature. The widely studied operators responsible for NSIs can be written as 



£ NS i = -2V2G F e f ^ c fc>fP L up) (f^Paf) , (6) 

f f 'c 

where sA are NSI parameters, a,p — e, /z,r, /, / = e,u } d, and C = L,R. Note that 
these operators are non-renormalizable and they are also not gauge invariant. Thus, using 
the NSI operators in Eq. ([6]), we find that the effective NSI parameters are 

2 

£ ^^—^^ (7) 
m x 

where my/ is the W boson mass and mx is the mass scale of the NSIs. If the new physics 
scale, i.e., the NSI scale, is of the order of 1(10) TeV, then one obtains e a/ 3 ~ 10~ 2 (10~ 4 ). 



In principle, NSIs can affect both (i) production and detection processes and (ii) prop- 
agation in matter, and one can have combination of both effects. In the following, we will 
first study production and detection NSIs, including the so-called zero-distance effect, and 
then matter NSIs. In addition, we will present mappings with NSIs and discuss approximate 
formulas for two neutrino flavors. 



A. Production and detection NSIs and the zero-distance effect 



In general, production and detection processes can be affected by NSIs, and for a realistic 
neutrino oscillation experiment, the neutrino states produced in a source and observed at a 
detector can be written as superpositions of pure orthonormal flavor eigenstates 
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-|65|: 



) = k>+ £ U\ u p) = (i + eWm) 

/3=e,/j.,T 



(8) 
(9) 



where the superscripts V and 'd' denote the source and the detector, respectively, and \u m ) 
is a neutrino mass eigenstate. Note that the states |z/*) and (vi\ are not orthonormal states 
due to the NSIs and that the matrices e s and e d are not necessarily the same matrix, since 
different physical processes take place at the source and the detector, which means that these 
matrices are arbitrary and non-unitary in general. If the production and detection processes 
are exactly the same process with the same participating fermions (e.g. /3-decay and inverse 
/3-decay), then the same matrix enters as e s = (z d )\ or on the form of matrix elements, 



£ a/3 = £ a/3 = ( £ ^a)* = ( £< j3a)* |66(| . For example, in the case of so-called non-unitarity effects 



which can be considered as a type of NSIs, see e.g. Ref. |67|) in the minimal unitarity 



violation model 
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-73], it holds that e s = Thus, it is important to keep in mind that 

these matrices are experiment- and pro cess- dependent quantities. 

In the case of production and detection NSIs, the neutrino transition probabilities are 



given by [65|, 
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In fact, an important feature of Eq. (jit) ft is that when a ^ (3, the first term is generally 
non-vanishing, which means that a neutrino flavor transition would already happen at the 
source before the oscillation process has taken place. This is known as the zero- distance 
effect 75|. It could be measured with a near detector close to the source. Note that Eq. (jTUj) 
is also usable to describe neutrino oscillations with a non-unitary mixing matrix, e.g. in the 
minimal unitarity violation model [68]. 



B. Matter NSIs 

In order to describe neutrino propagation in matter with NSIs, the simple effective mat- 
ter potential in Eq. ([3]) needs to be extended. Similar to standard matter effects, NSIs can 
affect the neutrino propagation by coherent forward scattering in Earth matter. The Earth 
matter effects are more or less involved depending on the specific terrestrial neutrino oscil- 
lation experiment. In other words, the Hamiltonian in Eq. fl3]) is replaced by an effective 
Hamiltonian, which governs the propagation of neutrino flavor states in matter with NSIs, 
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H = — [£/diag(r< m|, m\)U ] + diag(/4, 0, 0) + Ae m ] 



(12) 



where the matrix e m contains the (effective) NSI parameters e a p, which are defined as 



£a/3 - E af3 at 



f,P 



(13) 



with the parameters being entries of the Hermitian matrix e^ p and giving the strengths 
of the NSIs and the quantity Nf being the number density of a fermion of type /. Unlike 
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FIG. 1. Schematic pictures of standard matter effects (left picture) and matter non-standard 
neutrino interactions (right picture). 

e s and e d , e m is a Hermitian matrix describing NSIs in matter, where the superscript 'm' 
is used to distinguish matter NSIs from production and detection NSIs. Thus, for three 
neutrino flavors, we obtain 
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The "1" in the 1-1-element of the effective matter potential in Eq. (fT4"|) describes the interac- 
tion of electron neutrinos with left-handed electrons through the exchange of W bosons, i.e., 
the standard matter interactions, whereas the NSI parameters in the effective matter 
potential describe the matter NSIs. See Fig. [TJfor schematic pictures of standard and non- 
standard matter effects. Now, the effective Hamiltonian H in Eq. (fT2|) can be diagonalized 
using a unitary transformation, and one finds 



H = —Udiagim^m^ml)^ 



(15) 



where mj 



1,2,3) denote the effective mass-squared eigenvalues of neutrinos and U is 



the effective leptonic mixing matrix in matter. Explicit expressions for these quantities can 



be found in Ref. 



79| 



In the case of matter NSIs, for a constant matter density profile (which is close to reality 
for most long-baseline neutrino oscillation experiments), the neutrino transition probabilities 
are given by 



P{v a -)> z/g; L) 



(16) 



i=l 
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where L is the baseline length. Comparing Eq. (!T6|) with the formula for neutrino transition 
probabilities in vacuum, one arrives at the conclusion that there is no difference between 
the form of the neutrino transition probabilities in matter with NSIs and in vacuum if one 
replaces the effective parameters m 2 and U in Eq. ([TBI) by the vacuum parameters m| and 
U. The mappings between the effective parameters and the vacuum ones are sufficient to 
study new physics effects entering future long-baseline neutrino oscillation experiments (see 
Section Hi Bj) . The important point is the diagonalization of the effective Hamiltonian H and 
the derivation of the explicit expressions for the effective parameters. Now, using Eq. ([TBI) , 
we can express the neutrino oscillation probabilities in matter with NSIs (for a realistic 



experiment) as follows [79] 



P(u a ^u a ,L) = l-Aj2\U m U* aj \ 2 sm 2 (^^) , (17) 

i>j ^ ' 

/ \ / Am 2 L\ f Am 2 L\ 

P(u a up; L) = -4 J2 ^ {KfipAiU^) sin 2 I -^J- J - 8 J J[ sin f -^J- J (18) 

i>j ^ ' i>j ^ ' 

where (a, (3) run over (e,/i), (/x, r), and (r, e) and the quantity J7" is defined through the 
relation 



J 2 — |^ai| 2 |^j| 2 |f^aj| 2 |^/3i| 2 — - ^1 + | C^CKi | 2 1 C^SJ | 2 + |^ai| 2 |^/3i| 2 

-|f) m | 2 - \U^ 2 - \U a3 \ 2 - |^| 2 ) 2 . (19) 



C. Mappings with NSIs 



In Ref. 



791 ] . model-independent mappings for the effective masses with NSIs during prop- 



agation processes, i.e., in matter, were derived, which are given by 
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FIG. 2. Neutrino oscillation probabilities for the v e —> channel as functions of the neutrino 
energy E. We have set 5 = ir/2 and e er = 0.01, and all other matter NSI parameters are zero. 
Solid (black) curves are exact numerical results, dashed (red) curves are the approximative results, 
and dotted (blue) curves are results without NSIs. This figure has been adopted from Ref. [73 |. 

as well as model-independent mappings for the effective mixing matrix elements with NSIs, 
which are given by 



si 3 e 15 A(s 23 e e ^ + c 2 ze e r) 



U, 



(:3 



u, 



e2 



I- A 
A 



+ 



I - A 

+ c 23 £ eA1 - s 23 £ er ; 



^s 23 + -4 [c23£ MT + S23C23 (e^ - e TT ) - S23C23 (e^r + £* r )] 



(23) 
(24) 
(25) 



where a = Am^/Am^ and A = i/Am^. We observe from the explicit mappings (T20]) - 
(|25|) that the effective parameters can be totally different from the fundamental parameters, 
because of the dependence on A and the NSI parameters e a/ g. In addition, in Fig. [21 neutrino 
oscillation probabilities including the effects of NSIs for the electron neutrino-muon neutrino 
channel are plotted. It is found that the approximate mappings agree with the exact numer- 
ical results to an extremely good precision. Note that a singularity exists around 10 GeV 
due to the limitation of non-degenerate perturbation theory that has been used to derive 
the model-independent mappings. 



11 



D. Approximate formulas for two neutrino flavors 



The three-flavor neutrino evolution given in Eq. ( I14p is rather complicated and cum- 
bersome. Hence, in order to illuminate neutrino oscillations with NSIs, we investigate the 
oscillations using two flavors, e.g., v e and v T . In this case, we have the much simpler two- 
flavor neutrino evolution equation 



dL 
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2E 



.0 \ , / 1 + e ee e eT 
U I \U ] + A 

Am 2 



-er <-tt 




(26) 



where L is the neutrino propagation length that has replaced time in Eq. (1141) . Using 
Eq. (I26p . one can derive the two-flavor neutrino oscillation probability 

' Am 2 M V 
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where 8 m and Am M are the effective neutrino oscillation parameters when taking into 



account NSIs. These parameters are re 
8 and Am 2 and given by (see e.g. Ref. 



ated to the vacuum neutrino oscillation parameters 



W) 



(Am 2 7 ) 2 = [Am 2 cos(2#) - A (1 + e ee - e TT )] 2 + [Am 2 sin(2#) + 2Ae eT ] 2 , (28) 
Am 2 sin(2#) + 2Ae eT 
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For further discussion on approximate formulas for two neutrino flavors with NSIs, see 
Ref. £|. 



IV. THEORETICAL MODELS FOR NSIS 

In order to realize NSIs in a more fundamental framework with some underlying high- 
energy physics theory, it is generally desirable that it respects and encompasses the SM 
gauge group SU(3) x SU(2) x U(l). Note that the theoretical models presented in this 
section only represent a small selection, there exists many other models in the literature. In 
a toy model, including the SM and one heavy SU(2) singlet scalar field S, we can have the 



following interaction Lagrangian 82j 



A 5 nt = -KpL^^LpS , (30) 



12 



where the quantities \ aj 3 (a, (3 = e,fi,r) are elements of the coupling matrix A, L a is a 
doublet lepton field, and 02 is the second Pauli matrix. Now, integrating out the heavy field 
S, generates an anti-symmetric dimension-6 operator at tree level 83J, i.e., 

£ d=6,as = ( _ p ^ ^ pd%) (31) 

where ms is the mass of the heavy field S as well as Pl and Pr are left- and right-handed 
projection operators, respectively. Note that this is the only dimension-6 operator, which 
does not give rise to charged-lepton NSIs. 



A. Gauge symmetry invar iance 

At high-energy scales, where NSIs originate, there exists SU(2) x U(l) gauge symmetry 
invariance. In general, theories beyond the SM must respect gauge symmetry invariance, 



which implies strict constraints on possible models for NSIs (see e.g. Ref. 84|). Therefore, 
if there is a dimension-6 operator on the form 

then this operator will lead to NSI parameters such as e^. However, the above form for a 
dimension-6 operator must be a part of the more general form 

which involves four charged-lepton operators. Thus, we have severe constraints from experi- 
ments on processes like \i —> 3e, i.e., 

BR(n 3e) < 1(T 12 , 

which leads to the following upper bound on the above chosen NSI parameter 

e" < 10- 6 . 



Note that the above discussion is only valid for dimension-6 operators, and can be extended 
to operators with dimension equals to 8 or larger, but this will not be performed here. 
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(c) (d) 
FIG. 3. Tree-level diagrams with exchange of a heavy triplet Higgs in the type-II seesaw model. 



This figure has been adopted from Ref. 



B. Different see-saw models 



85[]. 



In a type-II seesaw model, the tree-level diagrams with exchange of a heavy triplet Higgs 
are given in Fig. [31 Integrating out the heavy triplet field (at tree level), we obtain the 



relations between the neutrino mass matrix and the NSI parameters as 



85| 



m, 



im, 



■)vf) ( m l) 



Op 



(32) 



a/3 " 8V2G F v% 

where v ~ 174 GeV is the vacuum expectation value of the SM Higgs field, wa is the mass 
of the Higgs triplet field, and is associated with the trilinear Higgs coupling. Now, using 
experimental constraints from lepton flavor violations and rare decays, we find upper bounds 
on the NSI parameters, which are presented in Table [Til From this table, we can observe 
that the NSI parameter e^ e has the largest upper bound. 

In addition, for m& = 1 TeV and varying m 1; we plot the upper bounds on some of the 
NSI parameters in the triplet seesaw model. The results are shown in Fig.Hl For a hierarchial 
mass spectrum (i.e., mi < 0.05 eV), all the NSI effects are suppressed, whereas for a nearly 
degenerate mass spectrum (i.e., mi > 0.1 eV), two NSI parameters can be sizable, which are 
e% and e™ 



C. The Zee Babu model 



In the Zee-Babu model 86M88| . we have the Lagrangian 

£ = £sm + f a pLl a Cia 2 L LI3 h + + g a pe^epk ++ - fxh~h~k ++ + h.c. + V H , (33) 
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TABLE II. Constraints on various e's from i — > ill, one-loop £ — > £j, and 
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FIG. 4. Upper bounds on NSI parameters in the type-II seesaw model. This figure has been 
adopted from Ref. 



where h + and k ++ are charged scalars that could be observed at the LHC and which lead to 
a two-loop diagram that generates small neutrino masses. The diagrams that are responsible 
for (a) NSIs of four charged leptons and (b) NSIs are presented in Fig. [51 

In Fig. El the allowed regions of the matter NSI parameters e MM and e TT in the Zee-Babu 
model are plotted. In the case of inverted neutrino mass hierarchy, the NSI parameters could 
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FIG. 5. Tree-level diagrams 
has been adopted from Ref. 
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FIG. 6. The allowed region of NSI parameters at la, 2a, and 3a C.L. in the Zee-Babu model. We 



take rrih = fn^ = /x = 10 TeV for the 
figure has been adopted from Ref. 



SSI- 



eft plot and m/i = rrik = fj, = 1 TeV for the right plot. This 



be in the range 10 4 4- 10 3 [89 1. 



V. PHENOMENOLOGY OF NSIS 

In this section, we will discuss the phenomenology of NSIs for atmospheric, accelerator, 
and reactor neutrino experiments as well as for accelerators, neutrino factory setups, and 
astrophysical settings such as solar and supernova neutrinos. As we will see, only two 
experimental collaborations have used their neutrino data to analyze NSIs, which are the 
Super-Kamiokande and MINOS collaborations. 
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A. Atmospheric neutrino experiments 



Neutrino oscillations with matter NSIs that are important for atmospheric neutrinos 
have been previously studied in the literature. For example, there are phenomenological 
studies that investigate the possibility to probe NSIs with atmospheric neutrino data only 

whereas there exist also more 



90 



92] and in combination with other neutrino data 



93 



theoretical investigations 



81 



991 ] . However, most importantly, there is an experimental study 



on m atter NSIs with atmospheric neutrino data from the Super-Kamiokande collaboration 



100J. In principle, atmospheric neutrinos are very sensitive to matter NSIs, since they travel 
over long distances inside the Earth before being detected 9^]. 

In order to analyze atmospheric neutrino data in the simplest way, we consider a two- flavor 
neutrino oscillation approximation with matter NSIs in the p^-Pt sector (cf. Section Till Dj) . 
since <H- v T oscillations are important for atmospheric neutrinos. In this case, the first 
row and first column in Eq. (|14p are cancelled, which effectively means that the NSI pa- 
rameters that couple to u e are set to zero, i.e., e ea = 0, where a = e, /i,r. In addition, 
the parameters Am^, 9±2, and 6^3 are not important, leading to a two-flavor approximation 



Am?,, 9 = # 23 , 



£ TP ,, and e^-r. In this 



that only includes the parameters Am = ^-w^ 31 , » — ^3, o MM , c MT - _ r/ ,. 
approximation, which has been named the two- flavor hybrid model, the two-flavor neutrino 
evolution equation reads 

. d 



dL 




1 
2E 



, \ . 1 + e 

U I I U ] + A 

Am 2 



fifj, ^fJbT 




(34) 



Using Eq. fl34l) . definin g e = £„ T and e' = e TT — and assuming that neutrinos have NSIs 



with c/-quarks only 90 



10l[ |, we obtain the two-flavor survival probability 91 



100] 



P{v„ ->■ 1/ ; L) = 1 - P{v„ ->• u T ; L) = 1 - sin 2 (26) sin 2 ( Am *} L R 



AE 



where the quantities and R are given by 

1 



sin 2 (26) = — [sin 2 (2#) + R 2 sin 2 (2£) + 2R sin(2#) sin(2£)] 



R = yjl + Rl + 2R [cos(2#) cos(2f) + sin(2fl) sin(2f)] 
with the two auxiliary parameters 



R = V2G F N f 



AE 



Am 2 

, 1 {2e 
£ — — arctan I — 



e' 2 

£ ' 2 + T 



(35) 

(36) 
(37) 

(38) 
(39) 
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FIG. 7. Allowed parameter regions for sin 2 (2^23) and AmJ 2 using the two-flavor hybrid model (solid 
curves) and standard two-flavor neutrino oscillations (dashed curves). The undisplayed parameters 



e and e' have been integrated out. This figure has been adopted from Ref. 



loot]. 



Now, using the two-flavor hybrid model together with atmospheric neutrino data from the 
Super-Kamiokande I (1996-2001) and II (2003-2005) experi ments , the Super-Kamiokande 
collaboration has obtained the following results at 90 % C.L. [100 ] 



M < l.i-nr 



and 



\e T r - < 4.9 • 1(T 



Amjj with and without 



and the allowed parameter regions for sin (2^23) and Am^ 
NSIs are shown in Fig. In principle, there are no significant differences between the 
allowed parameters regions with NSIs and the ones without NSIs. This analysis can be 
extended to a similar analysis with a three-flavor hybrid model also taking into account the 
NSI parameters e ee and e eT , which, however, leads to no significant changes for the allowed 
values of the parameter regions compared to the two-flavor hybrid model. It should be noted 
;hat the atmospheric neutrino data have no possibility to constrain the NSI parameter e ee 
95| and the other NSI paramet er e PT is related to both e ee and e TT via the expression 



e TT ~ 3|e eT | 2 /(l + 3e ee ) 2 [92 



95 



. In conclusion, the Super-Kamiokande collaboration 



has found no evidence for matter NSIs in its atmospheric neutrino data. 



B. Accelerator neutrino experiments 

Studies of previous, present, and future setups of accelerator neutrino oscillation experi- 
ments including matter NSIs have been thoroughly investigated in the literature, especially 

18 



setups with long-baselines belong to these studies. Such studies include s earc hes 



NSIs with the K2K experimen t 95(, the MINOS experiment 



80 



97 



99 



or matter 



102l - ll04j . and the 



OPERA experiment 105Hl07j ]. as well as sensitiv ity analyses of the NOi/A experiment 



and the T2K and T2KK experiments 



108 [ 



109h111| . The prospects for detecting N SIs at th e 



MiniBooNE experiment, which has a shorter baseline, has been inve stiga ted too 



112 



25 



98 



113 1 . 



114j . In addition, the 



There are also studies that are more general in character 
MINOS experiment has recently presented the results of a se arch for matter NSI in form of 
a poster at the "Neutrino 2012" conference in Kyoto, Japan 115 1 . 

Using three-flavor neutrino oscillations with matter NSIs for accelerator neutrinos, we will 
present the important NSI parameters and flavor transition probabilities for two experiments, 
which are (i) the MINOS experiment with baseline length L ~ 735 km (from Fermilab 
in Illinois, USA to Soudan mine in Minnesota, USA) and neutrino energy in the interval 
(1 -i- 6) GeV and (ii) the OPERA experiment with baseline length L ~ 732 km (from 
CERN in Geneva, Switzerland to LNGS in Gran Sasso, Italy) and average neutrino energy 
E ~ 17 GeV. Note that the baseline lengths of the two experiments are nearly the same, 
but there is a difference in the neutrino energy, which is about one order of magnitude. 

First, in the case of the MINOS experiment, the important NSI parameters are e eT and 
Eq. (TL4"j) ] and the interesting transition probability is the survival probability (or 
equivalently the disappearance probability), which to leading order is given by 



103] 



P(vp L) ~ 1 - sin^(2#2 3 ) sin 



-(Oil , .,;,.2 / ^ m 3l£ 



AE 



(40) 



where three-flavor effects due to Am^ and #13 have been neglected, and the effective para- 
meters are 



Am 31 = Am 3l£ > 



sin 2 (2M 



sin 2 (2fl 23 ) 



(41) 
(42) 



with 



As TT -\- cos(26>23) + sin 2 (2#23) 



(43) 



Note that in Eqs. (14T]) - (14"3~]) we have used the NSI parameter perturbation and the 

formulas sh ould hold if \e eT \ 2 A 2 L 2 /(AE 2 ) 1 or, in the case of the MINOS experiment, if 



\e eT \ < 5.8 



103]. Furthermore, it is possible to show the effective three- flavor mixing matrix 
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element U e3 is given by 

U e3 ~ U e3 + Ae eT cos(9 23 ), (44) 
which means t hat t here could be a degeneracy between the mixing angle 9\ 3 and the NSI 



parameter e t 



103[ | . Now, since the mixi ng; a ngle #13 has been measured 



116 



11% the 



MINOS experiment can put a limit on \e, 

MINOS experiment, it has recently been argued in the literature that it s 
to study the NSI parameter e^ T using the MINOS experiment too [97 



103]. In addition to the above discussion for the 

lould be possible 



99 



(assuming #23 = 45°), the survival probability becomes 99 ] 



1041 . In this case 



P{vn ->■ v» L) ~ 1 



sm 



Am| 1 



AE 



A 
2E 



(45) 



Now, using a model based on the survival probability in Eq. (I45p together with data from 
the MINOS experiment, the MINOS coll abor ation has obtained the following result for the 
matter NSI parameter e^ T at 90 % C.L. [1 15 ] 



-0.200 < e^ T < 0.070 



which means that MINOS has found no evidence for matter NSIs in its neutrino data, at 
least not a non-zero value for the matter NSI parameter e^ T . 

Second, in the case of the OPERA experiment, the important NSI parameter is £ Mr [see 
Eq. flHj) ] due to the relatively short baseline length and the interesting transiti on p robability 



is the appearance probability for oscillations of v„ into u T , which is given by 



107] 



P{y„ -»■ u T ; L) 



9 . „ s Am,, , A 



L 2 + 0(L 3 ) 



(46) 



where it has been assumed that the small mass-squared difference Am^ = 0. Thus, there 
exists a degeneracy between the fundamental neutrino oscillation parameters and the NSI 
parameter e^ T . Note that it has bee n sho wn that the OPERA experiment is not very sensitive 
to the NSI parameters e er and e TT [106 1 . 



C. Reactor neutrino experiments 



To my kn owle dge, NSIs in reactor neutrino experiments have only been discussed in 



Refs. |66|, 



1201 ] . Below, we will summarize these three works. 



20 



0.00 0.01 0.02 0.03 0.04 0.05 2 4 6 



FIG. 8. Mappings among #13, #13, and NSI parameters s a p. This figure has been adopted from 
Ref. m 



First, in Ref. 66[, a combined study on the performance of reactor and superbeam neu- 
trino experiments in the presence of NSIs is presented. Indeed, in this work, the authors 
argue that reactor and superbeam data can be used to establish the presence of NSIs. 

Second, in Ref. [74J, NSIs at reactor neutrino experiments only were studied. In Fig. [BJ 
mappings among the effective mixing angle #13 , the fundamental mixing angle #13, and 
the NSI parameters e a p are plotted. Such plots could be important for the analyses of, 
for example, the Daya Bay, Double Chooz, and RENO experiments. It was found that 
(i) #13 < 14°, which is larger than the former CHOOZ bound that is about 10° and also 
larger than the recently measured values of the mixing angle #13 and (ii) inspite of a very 
small #13 , a sizable effective mixing angle can be achieved due to mimicking effects [74] . In 
principle, this means that the measured value for the mixing angle #13 9° by the Daya Bay, 



Double Chooz, and RENO experiments 116l4ll9| could be a combination of the fundamental 



value for #13 (which should be smaller than the effective measured value) and effects of NSI 
parameters. 
Third, in Ref. 



120( | . NSIs at the Daya Bay experiment were studied. The authors of 
this work show that, under certain conditions, only three years of running of the Daya Bay 
experiment might be sufficient to provide a hint on NSIs. Thus, in future analyses of data 
from the Daya Bay experiment, it will be important to disentangle effects of NSI parameters 
on the mixing angle #13 as well as the large neutrino mass-squared difference Am^. In 
Fig. [9j the effects of NSIs on 13 and Am| 2 ~ Am^ are shown for Daya Bay after three 
years of running. 
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FIG. 9. Effects of NSIs on a nd Am^ for Daya Bay after three years of running. This figure 
has been adopted from Ref. [120]. 



D. Accelerators 



In Refs. 



121 



122], the authors have explored an alternative approach to study NSIs 



by bringing into play collider data. Indeed, LEP provides bounds on NSIs of the order of 
e < (10~ 2 ~r~ 10 -3 ) at y/s ~ 200 GeV, and if NSIs are contact interactions at LHC energies, 
LHC should have a sensitivity reach of NSIs of the order of e > 3 • 10 -3 at 14 TeV and with 
100 fb _1 of data. In addition, it should be mentioned that b ound s on the NSI parameters 
£-?f, e^, and (q = u,d and P — L,R) have been derived 1231 ] . using 1 ftT 1 data from 
the ATLAS experiment at LHC. 



E. Neutrino Factory 



Due to the large sensitivity at a future neutrino factory for small parameters such as NSI 
parameters, there exist a vast amount of inves tigations for different neutrino factory setups 



in connections to NSIs 



58 



62 



65. 



85 



89 



93 



124H131I ]. In this section, we will present 



most of these investigations and what their general conclusions for NSIs at a future neutrino 
factory are. Several investigations have discussed the sensitivity and discovery reach of a 



neutrino factory in probing NSIs 



3, 



124-127 



1311 ] . For example, it has been suggested that 



i) a 100 GeV neutrino factory could probe flavor-changing neutrino interactions of the order 



of |e| < 10" 4 at99% C.L. 



93| , (ii) there is an entanglement between the mixing angle #13 and 
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NSI para meters^ w hich would be solved in the best way by using the appearance channel 



130 



124 



125|, (iii) there are degeneracies between CP violation and NSI parameters 



13JJ , (iv) a neutrino factor y has excellent prospects in detecting NSIs originating from 



"new physics" at the TeV scale |126j ]. and (v) off-diagonal NSI parameters could be tested 
down to the order of 10~ 3 , whereas diagonal NSI parameter combinatio ns su ch as e ee — e TT 
and e w — e TT could only be tested down to 10 _1 and 1CT 2 , respectively [13l|. Furthermore, 
there are studies on the o ptimi zation of a neutrino factory with respect to NSI parameters, 



at a neutrino factory on NSIs 



especially for and e TT [128 1, as well as the impact of near detectors (with v T detection) 



1291 ] . In addition, it has been pointed out that the generation 



of matter NSIs might give rise to production and detection NSIs at a neutrino factory 
On the more theoretical side, some investigations on NSIs stemming from different 



models have been carried out. For example, NSIs from a triplet seesaw model |85j or the 



Zee-Babu model [89[ have been investigated. At a neutrino factory, NSIs from the triplet 
seesaw model could lead to quite significant signals of lepton flavor violating decays, whereas 
production NSIs from the Zee-Babu model might be at an observable level in the v e — > v T 
and/or — > v T channels. Finally, it has been pointed out that NSIs and non-unitarity 
might phenomenologically lead to very similar effects for a neutrino factory, although for 
completely different reasons 65 1. 



F. NSI effects on solar and supernova neutrino oscillations 



In addition to studies on NSIs with atmospheric and man-made sources of neutrinos, 
;here exi st some investigations of NSI effects on solar and supernova neutrino oscillations 



132 



1451 ] . First, we focus on NSIs with solar neutrinos, and then, we discuss NSIs with 



supernova neutrinos. 

In the case of solar neutrinos, an analysis of data from the Su per-K amiokande and SNO 



1331. Furthermore, it has 



experiments was performed to investigate the sensitivity of NSIs 
been suggested that the Borexino experiment can provi de signatures for NSIs |133| . and its 
data can also be used to place c onst raints on NSIs 1401 ] . The LENA proposal has also been 



illustrated as a probe for NSIs 



139] . In addition, non- universal flavor-conserving couplings 



1 This could be excluded due to the non-zero and quit e large v alue for the mixing angle #13 found by the 



Daya Bay, Double Chooz, and RENO collaborations 



116 11 



3- 



23 



with electrons, flavor-c hang i ng intera ctions, and NSIs in general can affect the phenomeno- 



logy of solar neutrinos 



132 



13j4l38|. 



In the case of supernova neutrinos, a three-flavor analysis for the possibility of probing 
NSIs, using neutrino s from a galactic supernova that propagate in the supernova envelope, 



has been performed 141] . Furthermore, the interp 



for supernova neutrinos has been investigated 



142 



ay 



between collective effects and NSIs 



to Ref. 
authors 



137 



144( | . Finally, a study on NSIs similar 



145]. 



for solar neutrinos has been presented for supernova neutrinos by the same 



In addition to NSIs with solar and supernova neutri nos, other studies on NSIs with 
astrophysics have been carried out. For example, in Ref. 1461 ] . the authors investigate pro- 
duction and detection NSIs of high-energy neutrinos at neutrino telescopes, using neutrino 
flux ratios. 



VI. PHENOMENOLOGICAL BOUNDS ON NSIS 

As discussed above, there are basically two neutrino experiments that have put direct 
bounds on NSI parameters - the Super-Kamiokande and MINOS experiments (see Sec- 
tions. IV Al and IV Bl) . However, there exist also some phenomenological works that have used 
different sets of data to find direct bounds on NSI parameters. Below, we will present direct 
bounds on both (i) matter NSIs and (ii) production and detection NSIs from such works. 
In addition, we will discuss NSIs for neutrino cross-sections. 



Direct bounds on matter NSIs 



In a work by Davidson et at 62j, bounds on matter NSI parameters using experiments 
with neutri nos and charged leptons have been derived for the realistic scenario with three 
flavors Id, ll27J, [l47] 



-0.9 < e ee < 0.75 \e efl \ < 3.8 ■ 10" 



kerl < 0.25 



\ 



V 



-0.05 < £ MjU < 0.08 |e MT | < 0.25 

M <0.4 



We observe that the bounds range from 10 4 to 1 for the different matter NSI parameters. 
Note that the bounds presented in this analysis are obtained using loop effects. However, it 
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turns ou t tha t bounds coming from loop effects are generally not applicable 14a ]. Therefore, 



in Ref. 



148 



Biggio, Blennow, and Fernandez- Martinez have performed a new analysis. 
The result of this analysis is that the model-independent bound for the NSI parameter 
e efl increases by a factor of 10 3 . It should b e no ted that one-loop effects on NSIs have 
also been studied in Ref. 
Refs. 



62 



135 



151 



1491. Now, in Ref. 



1501 ]. using bounds on NSI parameters from 



1521 ] . but disregarding the loop bound on the parameter e^, bounds on 



the effective matter NSI parameters have been estimated by Biggio, Blennow, and Fernandez- 
Martinez. The model-independent bounds on the matter NSI parameters (for the Earth) 
are given by 

\ 



< 2.5 e 



MM | 



< 0.21 

< 0.046 \e, 



< 1.7 

.,,r| < 0.21 

|e rT | < 9.0 

which, except for the matter NSI parameters and e MT , are larger than the too stringent 



V 



bounds found by Davidson et al. 



621 ]. and ranging between 10 2 and 10, i.e., they are one 



or two orders of magnitude larger than the previous bounds. 



B. Direct bounds on production and detection NSIs 



Finally, in Ref. 



150J, the authors have also derived model-independent bounds on pro- 



duction and detection NSIs. The most stringent bounds for charged-current-like NSIs for 
terrestrial experiments are the following 



e£| < 0.025 \e^\ < 0.030 |e£| < 0.030 
|ep < 0.025 \e^\ < 0.030 \e^ T \ < 0.030 
\e<i.t\ < 0.025 |e£f.| < 0.030 \e^\ < 0.030 



\ 



-TjU 



-ad 

'ee 




\e^\ < 0.025 



|e£| < 0.078 



/ 

< 0.04l\ 



\e^i\ < 0.013 



\ t TfJ,' 



\ett\ < 0.13 



/ 



0.013 
0.018 

where, whenever two values are presented, the upper value refers to left-handed NSIs and 
the lower one to right-handed NSIs. These bounds are basically of the order of magnitude 
lO" 2 . 
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C. NSIs for neutrino cross-sections 



In this section, NSIs for neutrino cross-sections will be investigated in detail. Neutrino 
NSIs with either electrons or first generation quarks can be constrained by low-energy scat- 
tering data. In general, one finds that bounds are stringent for muon neutrino interactions, 
loose for electron neutrino interactions, and in principle, do not exist for tau neutrino inter- 
actions. Note that in the present overview of the upper bou nds on the NSI parameters, the 
results from Biggio, Blennow, and Fernandez-Martinez 150| have not been included. 

The best measurement on electron neutrino-electron scatte ring comes from the LSND 
experiment, which found the cross-section of this process to be 153 1 

G zpTTlpE,, 



a(v e e -> ve) = (1.17 ±0.17)- 



(47) 



where m e is the electron mass and E v is the neutrino energy, which should be compared with 
the SM cross-section a(v e e — > ^e)| SM ~ 1.0967G ® F m e E v /'K. Including NSIs, the expression 



for this cross-section becomes 
2G 2 F m e E v 



a{y e e — >■ ve) 



Ti 



eL\2 



3^ 



_eB|2 



(48) 



where g e L ~ —0.2718 and g e R ~ 0.2326. Using the LSND data and the NSI cross-section for 
electron neutrino-electron scattering, we obtain 90 % CL. bounds on the NSI parameters 
(note that we are only considering one NSI parameter at a time) 



62] 



-0.07 < e!f < 0.11 , -l<£ff ? <0.5 



for flavor- conserving diagonal NSI parameters and 



|e^|<0.4, |^| < 0.7 



for flavor- changing NSI parameters. Considering both left- and right-handed diagonal NSIs, 
i.e., two NSI parameters simultaneously, a 90 % CL. region between two ellipses is obtained 



0.445 < (0.7282 + e^) 2 + -(0.2326 + e e *y < 0.725 

3 



which is shown in Fig. [TOJ Note that, in an update 154] . more restrictive allowed 
90 % CL. bounds for the NSI parameters and elf have been found 

-0.02 <e±< 0.09 - 0.11 <e e 3< 0.05. 
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FIG. 10. Bounds on flavor-conser ving NSIs of u e e scattering from the LSND experiment. This 



figure has been adopted from Ref. 



62] 



Then, we study electron neutrino-quark scattering. The CHARM collaboration 155| has 
measured the ratio between cross-sections and found 



R e 



o(v e N -> vX) + a(i> e N uX) 



cr(u e N ->■ eX) + a(u e N -» eX) 
Including NSIs, the quantities gie and can be expressed as 



(he) 2 + {QRef = 0.406 ±0.140. (49) 



(he) 2 = (si + efef + E l*«l 2 + (Sl + + El £ 



dL|2 
ae I ' 



a^e 



a^e 



(he) 2 = (9 U R + elff + E K«\ 2 + (9 d R + « 2 + E l^l" 



(50) 
(51) 



Using the CHARM data, we obtain 90 % C.L. bounds on the NSI parameters (only one NSI 
parameter at a time) 62] 



-1 < < e L < 0.3 , -0.3 < ef e < 0.3 , -0.4 < e£f < 0.7 , -0.6 < <f < 0.5 

for flavor-conserving NSIs and 

|4f|<0.5, q = u,d and P = L,R 

for flavor-changing NSIs. Again, a 90 % C.L. region considering several NSI parameters 
simultaneously can be obtained 

0.176 < (0.3493 + e^f + (-0.4269 + ef e f + (-0.1551 + e u e ff + (0.0775 + eff f < 0.636 . 
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Next, for muon neutrino-electron scattering, the CHARM II collaboration 156| has mea- 
sured g e v = -0.035±0.017 and g\ = -0.503±0.017, which translate into g\ = -0.269±0.017 
and g e R = 0.234 ±0.017. Using these values, one can compute 90 % C.L. bounds on the NSI 
parameters (only one NSI parameter at a time) 62] 



-0.025 < ef u < 0.03 



for flavor diagonal NSIs and 



K\ < o.i 



-0.027 < < 0.03 



L,R 



for flavor- changing NSIs. Furthermore, the NuTeV collaboration |157j has measured 
{dLfi) 2 = 0.3005 ±0.0014 and (gR^) 2 = 0.0310 ±0.0011 that appear in ratios of cross-sections 
for neutrino-nucleon (muon neutrino-quark) scattering processes. Using these values, we 
obtain (only one NSI parameter at a time) 62] 



-0.009 < < -0.003 



0.002 < ef u < 0.008 , 

dR 



-0.008 < < 0.003 , -0.008 < e a u n u < 0.015 



/'/' 



for flavor diagonal NSIs and 



< 0.05, q = u,d 



for flavor- changing NSIs. Similarly, the 90 % C.L. regions using two NSI parameters simul- 
taneously are presented in Fig. [TTJ 

Finally, we investigate NSIs for the e + e~ — > z/z/7 cross-section at LEP II. The 90 % C.L. bounds 
on flavor diagonal NSIs are given by 61] 



0.6<£^<0.4, -0.4 < 4? < 0.6, 



eR 

TT 



whereas for flavor-changing NSIs 



-ePi 

-a/31 



<0.4, P = L,R, a = r, p = e,fj,. 



Also, for the NSI parameters and e^, there exists an update |l54j of the 90 % C.L. bounds 

-0.51 < e± < 0.34 - 0.35 < e e * < 0.50 . 



In conclusion, using neutrino cross-section measurements for the determination of upper 
bounds on NSIs, the only non-zero parameters obtained are e u J^ and s'jz. Therefore, for 
future analyses, new data on neutrino cross-sections would be valuable. 
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FIG. 11. Bounds on flavor-conser ving NSIs of v^q scattering from the NuTeV experiment. This 



figure has been adopted from Ref. 



62] 



VII. SENSITIVITIES AND DISCOVERY REACH OF NSIS 



I am sure that more investigations about NSIs will be conducted in the future, both 
theoretically and experimentally. In addition, I believe that a future neutrino factorjo should 
be the most plausible experiment to find signatures of NSIs (see the detailed discussion in 
Section IV E|) . but perhaps also the LHC will shed some light. The power of a neutrino 
factory is that it may have a sensitivity and discovery reach for the NSI parameters. 



VIII. SUMMARY AND CONCLUSIONS 



In summary of this review, we have discussed NSIs as sub-leading effects to the standard 
paradigm for neutrino flavor transitions based on the phenomenon of neutrino oscillations. 
Especially, we have presented both (i) production and detection NSIs including the so-called 
zero-distance effect as well as (ii) matter NSIs. In the case of matter NSIs, we have given 
approximate analytical model-independent mappings between effective neutrino masses and 
leptonic mixing angles and the fundamental neutrino mass-squared differences and leptonic 
mixing parameters. In addition, we have studied approximate two-flavor formulas for neu- 
trino flavor transitions. Furthermore, we have presented some different theoretical models for 
2 Note that an interesting and important question is if a neutrino factory is going to be built or not, given 
a non-zero and quite large value for the mixing angle #13. 
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NSIs such as seesaw models and the Zee-Babu model, and investigated the phenomenology 
of NSIs in general. In particular, we have displayed the experimental results of upper bounds 
on NSIs from the analyses of the Super-Kamiokande (\e^ T \ < 1.1-10 -2 , \e TT — <4.9-lCT 2 ) 
and MINOS (—0.200 < £ MT < 0.070) collaborations, which both mean that no evidence for 
matter NSIs has been found. Moreover, we have indicated the sensitivities of NSIs for ac- 
celerators, a future neutrino factory, and the reactor neutrino experiment Daya Bay. For 
example, mimicking effects induced by NSIs could play a very important role in reactor 
neutrino experiments, especially for the mixing angle #13. In fact, the fundamental value 
for 9 13 could be smaller than the measured value due to an interplay with NSIs. Finally, we 
have presented phenomenological upper bounds on NSIs. In the case of matter NSIs, former 
results gave bounds ranging from 10~ 4 to 1, whereas latter results yielded bounds between 
10 -2 and 10. In the case of production and detection NSIs for terrestrial experiments, there 
are bounds of the order from 10 -6 to 0.1. In addition, we have given bounds for NSIs using 
data from neutrino cross-section measurements. Indeed, low-energy neutrino scattering ex- 
periments measuring neutrino cross-sections can be used to set bounds on NSI parameters. 
In conclusion, we have shortly discussed the future sensitivity and discovery reach of NSIs, 
which could be responsible for neutrino flavor transitions on a sub-leading level. Especially, 
the LHC and a future neutrino factory could open a new window towards determining the 
possible NSIs. 

ACKNOWLEDGMENTS 

I would like to thank Mattias Blennow and He Zhang for useful discussions and comments. 
This work was supported by the Swedish Research Council (Vetenskapsradet), contract no. 
621-2011-3985 (T.O.). 



[1] Y. Fukuda et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 81, 1562 (1998) 



arXiv:hep-ex/9807003 



[2] T. Ohlsson, Nature 485, 309 (2012). 



[3] L. Wolfenstein, Phys. Rev. D17, 2369 (1978) 



[4] S. Mikheyev and A. Y. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985). 



30 



[5] 



T. Schwetz, "Review on global fits," (2012), talk presented at the workshop "What's vl: 
Prom new experimental neutrino results to a deep understanding of theoretical physics and 
cosmology", Florence, Italy, June 25, 2012. 



D. Forero, M. Tortola, and J. Valle, (2012), arXiv:1205.4018 [hep-ph 
G. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo, et al, 
|Phys. Rev. D86, 013012 (2012)[|arXiv:1205.5254 [hep-phj[ 

Q. R. Ahmad et al. (SNO Collaboration), |Phys. Rev. Lett. 87, 071301 (2001)[ 



arXiv:nucl-ex/0106015, 



[9] 

[10] 

[11] 
[12] 

[13] 
[14] 

[15] 

[16] 
[17] 
[18] 

[19] 
[20] 

[21] 

[22] 



K. Eguchi et al. (KamLAND Collaboration), Phys. Rev. Lett. 90, 021802 (2003) 



arXiv:hep-ex/0212021 



M. Blennow, T. Ohlsson, and W. Winter, JHEP 06, 049 (2005), |aTXiv:hep-ph/0502147 
Y. Grossman and M. P. Worah, (1998), |arXiv:hep-ph /98075 1 1 

E. Lisi, A. Marrone, and D. Montanino, |Phys. Rev. Lett. 85, 1166 (2000) 



arXiv:hep-ph/0002053 



S. L. Adler, [Ph^sT Rev. D62, 117901 (2000)| |arXiv:hep-ph/00 05220. 



A. Gago, E. Santos, W. Teves, and R. Zukanovich Funchal, Phys. Rev. D63, 073001 (2001) 



arXiv:hep-ph/0009222, 



A. Gago, E. Santos, W. Teves, and R. Zukanovich Funchal, Phys. Rev. D63, 113013 (2001) 



arXiv:hep-ph/0010092 



T. Ohlsson, Phys. Lett. B502, 159 (2001), arXiv:hep-ph/0012272 



A. Gago, E. Santos, W. Teves, and R. Zukanovich Funchal, (2002), arXiv:hep-ph/0208166 



G. Fogli, E. Lisi, A. Marrone, and D. Montanino, Phys. Rev. D67, 093006 (2003) 



arXiv:hep-ph/0303064 



G. Barenboim and N. E. Mavromatos, |JHEP 050 1, 034 (2005)[ |arXiv:hep-ph/0404014 
D. Morgan, E. Winstanley, J. Brunner, and L. F. 



Thompson, 



Astropart. Phys. 25, 311 (2006), arXiv:astro-ph/0412618 



L. A. Anchordoqui, H. Goldberg, M. Gonzalez-Garcia, F. Halzen, D. Hooper, et al., 



Phys. Rev. D72, 065019 (2005), arXiv:hep-ph/0506168 



G. Barenboim, N. E. Mavromatos, S. Sarkar, 



and A. Waldron-Lauda, 



Nucl. Phys. B758, 90 (2006), arXiv:hep-ph/0603028 



31 



G. Fogli, E. Lisi, A. Marrone, D. Montanino, and A. Palazzo, 

|Phys. Rev. D76, 033006 (2007)[ |arXiv:0704.2568 [hep-ph][ 

J. Alexandre, K. Farakos, N. Mavromatos, and P. Pasipoularides, 



Phys. Rev. D77, 105001 (2008), arXiv:0712.1779 [hep-ph 



N. C. Ribeiro, H. Nunokawa, T. Kajita, S. Nakayama, P. Ko, et al, 



Phys. Rev. D77, 073007 (2008), arXiv:0712.4314 [hep-ph] 



N. E. Mavromatos, A. Meregaglia, A. Rubbia, A. Sakharov, and S. Sarkar, 



Phys. Rev. D77, 053014 (2008), arXiv:0801.0872 [hep-ph 



Y. Farzan, T. Schwetz, 



and A. Y. Smirnov, JHEP 0807, 067 (2008) 



arXiv: 0805. 2098 [hep-ph]. 



J. N. Bahcall, N. Cabibbo, and A. Yahil, Phys. Rev. Lett. 28, 316 (1972) 



V. D. Barger, W.-Y. Keung, and S. Pakvasa, Phys. Rev. D25, 907 (1982) 



J. Valle, Phys. Lett. B131, 87 (1983) 



V. D. Barger, J. Learned, S. Pakvasa, and T. J. Weiler, Phys. Rev. Lett. 82, 2640 (1999) 



arXiv:astro-ph/9810121 , 



S. Choubey and S. Goswami, Astropart. Phys. 14, 67 (2000), arXiv:hep-ph/9904257 



V. D. Barger, J. Learned, P. Lipari, M. Lusignoli, S. Pakvasa, et al., 



Phys. Lett. B462, 109 (1999), arXiv:hep-ph/9907421 



G. L. Fogli, E. Lisi, A. Marrone, and G. Scioscia, |Phys. Rev. D59, 117303~(1999) 



arXiv:hep-ph/9902267, 



S. Pakvasa, AIP Conf. Proc. 542, 99 (2000), arXiv:hep-ph/0004077 



S. Choubey, S. Goswami, 



and D. Majumdar, Phys. Lett. B484, 73 (2000) 



arXiv:hep-ph/0004193 



A. Bandyopadhyay, S. Choubey, and S. Goswami, Phys. Rev. D63, 113019 (2001) 



arXiv:hep-ph/0101273 



M. Lindner, T. Ohlsson, 



and W. Winter, Nucl. Phys. B607, 326 (2001) 



arXiv:hep-ph/0103170 



M. Lindner, T. Ohlsson, 



and W. Winter, Nucl. Phys. B622, 429 (2002) 



arXiv:astro-ph/0105309 



A. S. Joshipura, E. Masso, and S. Mohanty, Phys. Rev. D66, 113008 (2002) 



arXiv:hep-ph/0203181 



32 



J. F. Beacom and N. F. Bell, |Phys. Rev. D65, 113009 (2002)HarXiv:hep-ph/0204111[ 

A. Bandyopadhyay, S. Choubey, and S. Goswami, Phys. Lett. B555, 33 (2003) 



arXiv:hep-ph/0204173 



J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, 



Phys. Rev. Lett. 90, 181301 (2003), arXiv:hep-ph/0211305 



D. Indumathi, (2002), arXiv:hep-ph/0212038 



S. Ando, |Phys7 Lett. B570, 11 (2003)[ |arXiv:hep-ph/0307169 



G. Fogli, E. Lisi, A. Mirizzi, and D. Montanino, Phys. Rev. D70, 013001 (2004) 



arXiv:hep-ph/0401227 



S. Ando, [Phys7 Rev. D70, 033004 (2004)| |arXiv:hep-ph/0405200| 

S. Palomares-Ruiz, S. Pascoli, and T. Schwetz, j jHEP 0509, 048 (2005) 

|arXiv:hep-ph/0505216| 



D. Meloni and T. Ohlsson, |Phys. Rev. D75, 125017 (2007)| |arXiv:hep-ph /0612279, 
Gonzalez-Garcia 



M. 



and 



M. 



Maltoni, Phys. Lett. B663, 405 (2008) 



arXiv:0802.3699 [hep-ph 



M. Maltoni and W. Winter, JHEP 0807, 064 (2008) , arXiv:0803.2050 [hep-ph 



P. Mehta and W. Winter, JCAP 1103, 041 (2011)[ |arXiv:1101.2673 [hep-ph 



P. Baerwald, M. Bustamante, and W. Winter, (2012), arXiv: 1208.4600 [astro-ph.CO 



Y. Ashie et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 93, 101801 (2004) 



arXiv:hep-ex/0404034 



T. Araki et al. (KamLAND Collaboration), Phys. Rev. Lett. 94, 081801 (2005) 



arXiv:hep-ex/0406035. 



P. Adamson et al. (MINOS Collaboration), Phys. Rev. Lett. 101, 131802 (2008) 



arXiv:0806.2237 [hep-ex 



P. Adamson et al. (MINOS Collaboration), Phys. Rev. Lett. 106, 181801 (2011) 



arXiv:1103.0340 [hep-ex 



M. Blennow, T. Ohlsson, and W. Winter, |Eur. Phys. J. C49, 1023 (2007) 



arXiv:hep-ph/0508175 



S. Mikheyev and A. Y. Smirnov, Nuovo Cim. C9, 17 (1986) 



Y. Grossman, | Phys. Lett. B359, 141 (1995)] |arXiv:hep-ph/9507344[ 



Z. Berezhiani and A. Rossi, Phys. Lett. B535, 207 (2002), arXiv:hep-ph/0111137 



33 



S. Davidson, C. Peha-Garay, N. Rius, and A. Santamaria, JHEP 03, Oil (2003), 
|arXiv:hep-ph/0302093 [ 

M. Gonzalez-Garcia, Y. Grossman, A. Gusso, and Y. Nir, |Phys. Rev. D64, 096006 (2001)[ 



arXiv:hep-ph/0105159 



S. M. Bilenky and C. Giunti, [Phys. Lett. B300, 137 (1993)} |arXiv:hep^ph /9211269 
D. Meloni, T. Ohlsson, W. Winter, and H. Zhang, JHEP 1004, 041 (2010) 



arXiv:0912.2735 [hep-ph 



J. Kopp, M. Lindner, T. Ota, and J. Sato, Phys. Rev. D77, 013007 (2008) 



arXiv:0708.0152 [hep-ph 



Z.-z. Xing, Int. J. Mod. Phys. A23, 4255 (2008), arXiv:0810.1421 [hep-ph 



S. Antusch, C. Biggio, E. Fernandez-Martinez, M. Gavela, and J. Lopez-Pavon, 
JHEP 0610, 084 (2006)1 |arXiv:hep-ph/0607020| 

E. Fernandez-Martinez, M. Gavela, J. Lopez-Pavon, and O. Yasuda, 
Phys. Lett. B649, 427 (2007)] |arXiv:hep-p£70703098 



S. Goswami and T. Ota, Phys. Rev. D78, 033012 (2008), arXiv:0802.1434 [hep-ph 



Z.-z. Xing and S. Zhou, Phys. Lett. B666, 166 (2008) , arXiv:0804.3512 [hep-ph 



S. Luo, |Phys. Rev. D78, 016 006 (2008)] [arXlv:0804.4897 [hep-phj 



G. Altarelli and D. Meloni, Nucl. Phys. B809, 158 (2009), arXiv:0809.1041 [hep-ph 



T. Ohlsson and H. Zhang, |Phys. Lett. B671, 99 (2009) , arXiv:0809.4835 [hep-ph 



P. Langacker and D. London, [Phys. Rev. D38, 907 (1988) 



J. Valle, Phys. Lett. B199, 432 (1987) 



M. Guzzo, A. Masiero, and S. Petcov, [Phys. Lett. B260, 154 (1991) 
E. Roulet, |Phys. Rev. D44, 935 (1991) 



D. Meloni, T. Ohlsson, and H. Zhang, JHEP 04, 033 (2009)[ [arXiv: 090 1.1784 [hep-ph] 



N. Kitazawa, H. Sugiyama, and O. Yasuda, (2006), |arXiv:hep-ph/0606013[ 



M. Blennow and T. Ohlsson, Phys. Rev. D78, 093002 (2008) , arXiv: 0805. 2301 [hep-ph 



M. S. Bilenky and A. Santamaria, |Nucl. Phys. B420, 47 (1994)] |arXiv:hep-ph/93l'0302| 



S. Antusch, J. P. Baumann, and E. Fernandez-Martinez, Nucl. Phys. B810, 369 (2009) 



arXiv:0807.1003 [hep-ph 



[84] D. Hernandez, (2009), arXiv:0911.4764 [hep-ph 



34 



[85] M. Malinsky, T. Ohlsson, and H. Zhang, |Phys. Rev. D79, 011301(R) (2009) 



arXiv:0811.3346 [hep-ph 



[86 
[87; 
[88 
[89 

[90; 

[91 

[92 

[93 
[94 

[95 
[96 

[97 



[99 
[100 
[101 
[102 



A. Zee, |Phys. Lett. B161, 141 (1985) 



A. Zee, Nucl. Phys. B264, 99 (1986) 



K. S. Babu, Phys. Lett. B203, 132 (1988) 



T. Ohlsson, T. Schwetz, and H. Zhang, |Phys. Lett. B681, 269 ( 2009) 

|arXiv:0909.0455 [hep-ph]} 



N. Fornengo, M. Maltoni, R. Tomas, and J. Valle, |Phys. Rev. D65, 013010 ( 2002) 



arXiv:hep-ph/0108043 



M. Gonzalez-Garcia and M. Maltoni, Phys. Rev. D70, 033010 (2004) 



arXiv:hep-ph/0404085 



A. Friedland, C. Lunardini, and M. Maltoni, Phys. Rev. D70, 111301 (2004) 



arXiv:hep-ph /0408264, 



P. Huber and J. Valle, Phys. Lett. B523, 151 (2001), arXiv:hep-ph/0108193 



M. Guzzo, P. de Holanda, M. Maltoni, H. Nunokawa, M. Tortola, et al, 



Nucl. Phys. B629, 479 (2002), arXiv:hep-ph/0112310 



A. Friedland and C. Lunardini, IPhys. Rev. D72, 053009 (2005), arXiv:hep-ph/0506143 



M. Gonzalez-Garcia, M. Maltoni, and J. Salvado, JHEP 1105, 075 (2011) 



arXiv:1103.4365 [hep-ph 



J. Kopp, P. A. Machado, and S. J. Parke, | Phys. Rev. D82, 113002 (2010)[ 



arXiv: 1009.0014 [hep-ph 



F. Escrihuela, M. Tortola, J. Valle, and O. Miranda, Phys. Rev. D83, 093002 (2011) 



arXiv:1103.1366 [hep-ph 



W. A. Mann, D. Cherdack, W. Musial, and T. Kafka, Phys. Rev. D82, 113010 (2010) 



arXiv: 10 06.5720 [hep-ph]] 



G. Mitsuka et al. (Super-Kamiokande Collaboration), Phys. Rev. D84, 113008 (2011) 



arXiv:1109.1889 [hep-ex 



M. Gonzalez-Garcia, M. Guzzo, P. Krastev, H. Nunokawa, O. Peres, et al., 



Phys. Rev. Lett. 82, 3202 (1999), arXiv:hep-ph/9809531 



A. Friedland and C. Lunardini, Phys. Rev. D74, 033012 (2006), arXiv:hep-ph/0606101 , 



35 



[103] M. Blennow, T. Ohlsson, 
arXiv:hep-ph/0702059 , 



and J. Skrotzki, |Phys. Lett. B660, 522 (2008) 



[104 
[105 
[106 

[107; 

[108 
[109 
[110 
[111 
[112 

[113 
[114 

[115 



[116 

[117; 

[118 
[119 

[120 

[121 
[122 



Z. Isvan (MINOS Collaboration), (2011), arXiv:1110.1900 [hep-ex 



T. Ota and J. Sato, Phys. Lett. B545, 367 (2002), arXiv:hep-ph/0202145 
A. Esteban-Pretel, J. Valle, and P. Huber, |Phys. Lett. B668, 197 (2008) 



arXiv: 0803. 1790 [hep-ph 



M. Blennow, D. Meloni, T. Ohlsson, F. Terranova, and M. Westerberg, 



Eur. Phys. J. C56, 529 (2008) , arXiv:0804.2744 [hep-ph 



A. Friedland and I. M. Shoemaker, (2012), arXiv:1207.6642 [hep-ph 



H. Oki and O. Yasuda, Phys. Rev. D82, 073009 (2010), arXiv: 1003.5554 [hep-ph 



O. Yasuda, Nucl. Phys. Proc. Suppl. 217, 220 (2011), arXiv:1011.6440 [hep-ph 



R. Adhikari, S. Chakraborty, A. Dasgupta, and S. Roy, (2012), arXiv: 1201 .3047 [hep-ph 
L. M. Johnson, R. S. Williams, L. F. Spencer, and B. J. DeWilde, (2006), 



arXiv:hep-ph/0610011 



E. Akhmedov and T. Schwetz, JHE P 1010, 115 (2010)| |arXiv: 1007.4171 [hep-ph 
T. Ota, J. Sato, 



and N.-a. Yamashita, Phys. Rev. D65, 093015 (2002) 



arXiv:hep-ph/0112329 



J. A. B. Coelho (MINOS Collaboration), "Search for effects of exotic models in MINOS," 
(2012), poster presented at the XXVth International Conference on Neutrino Physics & 
Astrophysics (Neutrino 2012), Kyoto, Japan, June 3-9, 2012. 



F. An et al. (Daya Bay Collaboration), Phys. Rev. Lett. 108, 171803 (2012) 



arXiv:1203.1669 [hep-ex 



Y. Abe et al. (Double Chooz Collaboration), Phys. Rev. Lett. 108, 131801 (2012) 



arXiv:1112.6353 [hep-ex 



Y. Abe et al. (Double Chooz Collaboration), (2012), arXiv: 1207.6632 [hep-ex 



J. Ahn et al. (RENO collaboration), Phys. Rev. Lett. 108, 191802 (2012) 



arXiv: 1204.0626 [hep-ex 



R. Leitner, M. Malinsky, B. Roskovec, and H. Zhang, JHEP 1112, 001 (2011) 



arXiv:1105.5580 [hep-ph 



S. Davidson and V. Sanz, Phys. Rev. D84, 113011 (2011), arXiv:1108.5320 [hep-ph 



S. Davidson and V. Sanz, (2011), arXiv:1110.1558 [hep-ph] 



36 



[123 

[124 

[125 
[126 

[127; 

[128 
[129 
[130 

[131 

[132 

[133 

[134 
[135 

[136 
[137 
[138 
[139 

[140 
[141 



A. Friedland, M. L. Graesser, I. M. Shoemaker, andL. Vecchi, |Phys. Lett. B714, 267 (2012) 
|arXiv: 1111.5331 [hep-ph]} 

P. Huber, T. Schwetz, and J. Valle, |Phys. Rev. Lett. 88, 101804 (2002)[ 



arXiv:hep-ph/0111224 



P. Huber, T. Schwetz, and J. Valle, Phys. Rev. D66, 013006 (2002) , arXiv:hep-ph/0202048 



J. Kopp, M. Lindner, and T. Ota,|Phys. Rev. D76, 013001 (2007)] |ar35v-hep-ph/0702269 



N. Ribeiro, H. Minakata, H. Nunokawa, S. Uchinami, and R. Zukanovich-Funchal, 



JHEP 0712, 002 (2007), arXiv:0709.1980 [hep-ph 



J. Kopp, T. Ota, and W. Winter, Phys. Rev. D78, 053007 (2008) , arXiv:0804.2261 [hep-ph 



J. Tang and W. Winter, |Phys. Rev. D80, 053001 (2009)[ |arXiv:0903.3039 [hep-ph]} 

A. Gago, H. Minakata, H. Nunokawa, S. Uchinami, and R. Zukanovich Funchal, 

JHEP 1001, 049 (2010)[ |arXiv:0904.3360 [hep-ph]] 



P. Coloma, A. Donini, J. Lopez-Pavon, and H. Minakata, JHEP 1108, 036 (2011) 



arXiv:1105.5936 [hep-ph 



S. Bergmann, M. Guzzo, P. de Holanda, P. Krastev, 



and H. Nunokawa, 



Phys. Rev. D62, 073001 (2000) [ |arXiv:hep-ph /0004049 



Z. Berezhiani, R. Raghavan, 



and A. Rossi, Nucl. Phys. B638, 62 (2002) 



arXiv:hep-ph/0111138 



M. Guzzo, P. de Holanda, and O. Peres, |Phys. Lett. B591, 1 (2004)[|arXiv:hep-ph /0403134 
A. Bolahos, O. Miranda, A. Palazzo, M. Tortola, 



and 



Valle, 



Phys. Rev. D79, 113012 (2009), arXiv:0812.4417 [hep-ph] 



A. Palazzo and J. Valle, Phys. Rev. D80, 091301 (2009), arXiv:0909.1535 [hep-ph 



C. Das and J. Pulido, Phys. Rev. D83, 053009 (2011), arXiv:1007.2167 [hep-ph 



A. Palazzo, Phys. Rev. D83, 101701 (2011), arXiv: 1101. 3875 [hep-ph 



E. Garces, O. Miranda, M. Tortola, and J. Valle, Phys. Rev. D85, 073006 (2012) 



arXiv:1112.3633 [hep-ph 



S. K. Agarwalla, F. Lombardi, and T. Takeuchi, (2012), [axXivT 1 207.349 2 [hep-ph 
A. Esteban-Pretel, R. Tomas, and J. Valle, |Phys. Rev. D76, 053001" (2007) 



arXiv:0704.0032 [hep-ph 



[142] M. Blennow, A. Mirizzi, and P. D. Serpico, Phys. Rev. D78, 113004 (2008) 



arXiv:0810.2297 [hep-ph 



37 



[143] 

[144] 

[145] 
[146] 
[147] 

[148] 

[149] 

[150] 

[151] 

[152] 

[153] 

[154] 
[155] 
[156] 
[157] 



and J. Valle, |Phys. Rev. D81, 063003 "(2010) 



A. Esteban-Pretel, R. Tomas, 
|arXiv:0909.2196 [hep-ph]] 

B. Dasgupta, G. G. Raffelt, and I. Tamborra, |Phys. Rev. D81, 073004 (2010) [ 



arXiv:1001.5396 [hep-ph 



C. Das and J. Pulido, (2011), arXiv: 1111.6939 [hep-ph 



M. Blennow and D. Meloni, |Phys. Rev. D80, 065009 (2009) [ |arXiv:0901.2110 [hep-ph] 
J. Abdallah et al. (DELPHI Collaboration), |Eur. Phys. J. C38, 395 (2005) 



arXiv:hep-ex/0406019. 



C. Biggio, M. Blennow, and E. Fernandez-Martinez, JHEP 03, 139 (2009) 



arXiv:0902.0607 [hep-ph 



B. Bellazzini, Y. Grossman, I. Nachshon, and P. Paradisi, JHEP 1106, 104 (2011) 



arXiv: 1012.3759 [hep-ph 



C. Biggio, M. Blennow, and E. Fernandez-Martinez, JHE P 08, 090 (2 009) , 



arXiv:0907.0097 [hep-ph 



J. Barranco, O. Miranda, C. Moura, and J. Valle, |Phys. Rev. D73, 113001~(2006) 



arXiv:hep-ph/0512195 



J. Barranco, O. Miranda, C. Moura, and J. Valle, Phys. Rev. D77, 093014 (2008) 



arXiv:0711.0698 [hep-ph 



L. B. Auerbach et al. (LSND Collaboration), Phys. Rev. D63, 112001 (2001) 



arXiv:hep-ex/0101039 



D. Forero and M. Guzzo, [Phys. Rev. D84, 013002 (2011) 



J. Dorenbosch et al. (CHARM Collaboration), |Phys. Lett. B180, 303 (1 986). 



P. Vilain et al. (CHARM-II Collaboration), Phys. Lett. B335, 246 (1994) 



G. P. Zeller et al. (NuTeV Collaboration), Phys. Rev. Lett. 88, 091802 (2002) 



arXiv:hep-ex/0110059 



38 



